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The giant flare of December 27, 2004 from SGR 1806-20 represents one of 
the most extraordinary events captured in over three decades of monitoring the 

q ! 7-ray sky. One measure of the intensity of the main peak is its effect on X- and 

\ 7-ray instruments. RHESSI, an instrument designed to study the brightest solar 

flares, was completely saturated for ~0.5 s following the start of the main peak. 
A fortuitous alignment of SGR 1806-20 near the Sun at the time of the giant flare, 
however, allowed RHESSI a unique view of the giant flare event, including the 
precursor, the main peak decay, and the pulsed tail. Since RHESSI was saturated 

.{ph ! during the main peak, we augment these observations with Wind and RHESSI 

particle detector data in order to reconstruct the main peak as well. Here we 
present detailed spectral analysis and evolution of the giant flare. We report 
the novel detection of a relatively soft fast peak just milliseconds before the main 
peak, whose timescale and sizescale indicate a magnetospheric origin. We present 
the novel detection of emission extending up to 17MeV immediately following 
the main peak, perhaps revealing a highly-extended corona driven by the hyper- 
Eddington luminosities. The spectral evolution and pulse evolution during the 
tail are presented, demonstrating significant magnetospheric twist and evolution 
during this phase. Blackbody radii are derived for every stage of the flare, which 
show remarkable agreement despite the range of luminosities and temperatures 
covered. Finally, we place significant upper limits on afterglow emission in the 
hundreds of seconds following the giant flare. 
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Introduction 



The soft gamma repeater SGR 1806-20 was discovered in 1979 (lLaros et al.lll986l ). and 
has been studied intensively over the intervening two decades at X-ray, 7-ray, infrared, 
and radio wavelengths. It has emitted over 450 soft 7-ray bursts, mostly of short dura- 
tion, during sporadic active periods, and has been found to be a quiescent, variable X -ray 
source as well, emitting up to ~ 150 keV jMereghetti et al.ll2005cl : iMolkov et al.ll2005l ). In- 
deed, X-ray observations of its periodic, quiescent c omponent have provided some of the 
best evide nce for a magnetar-strength m agne tic field (|Kouveliotou et al.lll998l ). as first pro- 
posed by iDuncan fc Thompson I (119921 ) and iPaczyhski I (Il992l ). The infrared counterpart 
to SGR 1806-20 is a faint, highly obscured source, in keeping with its location towards the 
Galactic center (IKosugi et al.ll2005l ; llsrael et al.ll2005l ). Presumably, it is a lone neutron star, 
whose infrared intensity varies roughly in concert with bursting activity and its quiescent 
X-ray flux. There have been numerous attempts to determine t he distance to SGR 1 806- 
20 by va rious methods, leading to estimates from 6.4 - 9.8 kpc (ICameron et al. 1 120051) to 



15.1 kpc flMcClure-Griffiths fc Gaenslerll2005l : Eikenberrv et al. 1 12004 ; ICorbel fc Eikenberrv 



20041 ; ICorbel et al. 1119971 ). In this paper, we will quote all energies and luminosities in terms 
of dio = (d/lOkpc). Like SGR 0525-66 and SGR 1900+14, SGR 1806-20 has emitted a long 
duration, hard spectrum gi ant flare, whose flux at Earth greatly exceeded that of any other 



known cosmic X-ray source (IHurley et al. 1120051 ; iMazets et al. 1120051 ; iMereghetti et al. 



1999 



Palmer et al.ll2005l ). SGRs are not detectable quiescent radio emitters (ILorimer fc Xilouris 



2005a; 



2000) , but giant flares creat e transient radio nebulae which are observable for weeks (IFrail et al 



Gaensler et all 120051 ). 



In the ma gnetar model, magnetic dissipatio n, rather than rotation, provides the main 
energy source ([Thompson fc Duncanlll995l . Il996l ). Steady dissipation heats the neutron star 
surface and powers the quiescent X-ray emission. Localized crustal cracking causes short- 
duration, soft spectrum, and relatively weak bursts during active periods. Major crustal 
reconfigurations are thought to be responsible for the rarer long-duration, hard spectrum, 
intense giant flares. 

In 2004, SGR 1806-20 underwent a period of intense activity. The rate of small burst s 
peaked around mid- year, in conjunction with the quiescent X-ray flux (IWoods et al. 1 120061 ). 
The spindown rate, as evidenced by the frequency derivative, decreased. This activity cul- 
minated in the giant flare of December 27, 2004. Unlike the case of the giant flare from 
SGR 1900+14, however, there was no sudden change to the spin frequency. X-ray observa- 
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tions carried out several months later, however, revealed a slower spin-down rate, a smaller 
pulsed fraction for the quiescent emission, a differen t pulse profile, a softer spectrum, and a 
decreased flux (ITiengo et al.ll2005t iRea et al. M2005I ) . In the magnetar model, these changes 
are attributed to a major reconfiguration of the neutron star's magnetic field. 

In this paper, we present a detailed analysis of the Ramaty High Energy Solar Spectro- 
scopic Imager (RHESSI) data on the giant flare, concentrating on the time-resolved energy 
spectra of its various phases from 3 keV to 17 MeV. By virtue of the high time and energy res- 
olution and broad spectral coverage of the measurements, we believe that these data present 
the most complete spectral picture of this, or any other giant flare. 



2. RHESSI Spectrometer 

RHESSI is an array of nine coaxial germanium detectors, designed to perform de- 
tailed spectrosco pic imaging of X-ray and 7-ray emission (3keV - 17 MeV) from solar flares 



(jLin et al.ll2002l ). Spectral resolution ranges from ~lkeV FWHM in the hard X-ray range, 



up to several keV in the MeV range. RHESSI imaging is performed by two arrays of opaqu e 



1-D grids, separated by 1.55 m, and co-aligned with the nine detectors (IZehnder et al.ll2003l ). 
As the RHESSI spacecraft rotates (4.07 s period, axis aligned with the Sun) these grids mod- 
ulate the count rat e in the detectors, a llowing imaging through rotational modulation colli- 



mator techniques (IHurford et al. 1 120021 ) . Thus, RHESSI has high angular resolution (2.3") 
in the 1° field of view of its optics. However, the detectors themselves are unshielded, and 
are able to view transient sources from the whole sky. In addition, RHESSI sends down the 
energy and timing information for each photon, allowing detailed timing measurements. 

RHESSI observed both the precursor to and the giant flare from SGR 1806-20 in their 
entirety, starting at 21:28:03.44 UT and 21:30:26.64 UT 2004-12-27 at the spacecraft respec- 
tively. At the time of this event, SGR 1806-20 was located 5° from the Sun, just outside the 
primary imaging field of view of the RHESSI instrument. At this angle from the solar direc- 
tion, the shadow pattern of one front grid falls on the aligned bottom grid of a neighboring 
detector once per rotation, a fortuitous alignment that allows us to get a ^-second "snap- 
shot" of the direct spectrum, down to 3 keV, twice per RHESSI rotation period (4.07 s). 
During the main peak of the flare the RHESSI spectroscopy detectors were saturated for 
~0.5s after the initial rise of the main peak, but observed the decay of the main peak and 
the 400-s long pulsed tail. 



Fig. 1 shows the 20-100 keV lightcurve of the SGR 1806-20 giant flare, from just before 
the precursor to after the end of the pulsed tail. In this plot and throughout this paper, we are 
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quoting times t26 in seconds relative to 21:30:26 UT 2004-12-27. (So the main peak starts at 
^26 = 0.65 s.) Since SGR 1806-20 was located just outside of the main RHESSI field-of-view, 
this lightcurve is dominated by photons which have scattered into the detectors from other 
parts of the spacecraft or have been Compton reflected from Earth's atmosphere. Therefore, 
we do not attempt to use these events for spectral analysis below ~0.2MeV. All of our hard 
X-ray spectra are derived using the snapshot data (with the exception of the fast peak), 
where RHESSI has a direct view of the flare. At higher energies, the RHESSI grids are 
more transparent and Earth reflection is negligible, allowing all photon events to be used for 
spectral analysis. 



For the snapshot spectra we have used the on-axis RHESSI response matrices (ISmith et al 



20021 ). which produce acceptable spectral fits to the snapshot data. For the RHESSI snap- 
shot spectra analyzed in this paper we assume an absorption of Nh = 6 .69 x 10 22 cm~ 2 , 



which was measured for this source prior and su bsequent to the giant flare (IMurakami et al. 
1994! : iMereghetti et al.ll2005bl : iRea et al-lliooi ). 



3. RHESSI & Wind Charged Particle Detectors 



During the intense main peak all X- and 7-ray instruments experienced some degree of 
saturation, making reliable reconstruction of the time history and energy spectrum difficult or 
impossible. Many particle detectors, however, are small, thin silicon detectors with very low 
effective areas for X- and 7-ray interactions. Most of these detectors are usually impervious 
to 7-ray photons; however, due to the brightness of the main peak a number of particle 
instruments reg istered strong signals without saturating, allowing detailed reco nstruction of 
the main peak (IHurley et al.ll2005l ; iTerasawa et al. 1 120051 ; ISchwartz et al. 1 120051 1 . 



The RHESSI particle detector (PD), which is used for detecting SAA passages, is a 
small silicon detector of area 0.25 cm 2 , 960 microns thick (ISmith et al. 1120021 ). The RHESSI 
PD is instrumented to measure simple counting rates for all interactions above two threshold 
levels, ~50keV and ~620keV, with 0.125-s time resolution. The PD, which normally does 
not register even the brightest solar flares, measured a strong burst of counts from the main 
peak (Fig. 2), but still well below saturation levels. Therefore, the PD gives us a very good 
measure of the incident count rates above these two thresholds. However, the limited PD 
data do not allow us to strongly constrain the shape of the 7-ray spectrum. 

Therefore, we al so analyzed the observations from the Wind 3D Plasma & Energetic 
Particle Experiment (jLin et al.lll995l ). Wind has six double-ended Solid State Telescopes 



(SSTs), five with two back-to-back 1.5 cm 2 , 300 micron thick silicon detectors (called O and 
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F, with 9 and 7 PHA channels, respectively), and one SST with a third, 1.5 cm 2 , 500 micron 
thick detector (T) in between. The multi-channel analyzers covered the 20keV - HMeV 
range with various time resolutions between 12 and 96 s. All of the Wind detectors see a 
strong signal from the main peak, and they were all used in our analysis. The Wind SST 
data is saturated, but in a fairly benign way. The detectors and shaping electronics have 
a fast response (0.5 /is shapers), and are well below saturation level. However, the analog- 
to-digital converters (ADC) for the telescopes are shared, and exceeded their maximum 
throughput during the main peak. When the ADC is busy processing an event, incoming 
events are thrown away, but without pileup (D. Curtis, private communication). Therefore, 
the Wind SSTs accurately sampled the spectral shape of the main peak, but not the overall 
normalization. 



Between the Wind SSTs and the RHESSI PD we are able to reconstruct the input 
spectrum and the overall normalization, respectively. While our work focuses on the spectral 
analysis, other particle detector observations have been able to reconstruct the main pea k 



lightcurve with much higher time resolution (ITerasawa et al. 



2005 



The summary of these spectral results was presented in lHurley et al. 
more details of the analysis method. 



Schwartz et al.ll2005l ). 



(120051 ) . Here we present 



4. Stages of the Giant Flare 

For the purpose of our analysis, the giant flare from SGR 1806-20 can be divided into 
six separate stages, illustrated in Fig. 1. There is a precursor flare (i) 142 seconds prior 
to the main peak, followed by a quiescent preflare period (ii). Immediately prior to the 
main peak, a fast peak (iii) occurs which lasts merely 2.5 ms. The main peak itself (iv) 
lasts ~0.5s. Subsequent to the main peak, there is a brief ~60s decay period dominated 
by strong nonthermal emission (v), followed by the characteristic pulsed tail (vi) lasting 
400 seconds. Finally, there is the postflare period with the potential for afterglow emission 
(vii). While RHESSI was saturated ~0.5s during the main peak, it had an excellent view 
of SGR 1806-20 during the rest of these stages. Here we present a detailed analysis of each 
stage of this spectacular event. 



4.1. Precursor 



RHESSI observed a precursor burst during 21:28:03.44-21:28:04.49 UT, 142 s before the 
main peak of the giant flare, with a peak count rate in the spectroscopy detectors of ~30,000 



- 6- 



cnts -1 . The relatively long duration (Is) and nearly flat lightcurve of this burst (Fig. 3) 
distinguish it from more common SGR bursts with typical durations of ~0.1s. We see a 
rise time for the precursor of 27 ms, and a fall time of 110 ms. RHESSI caught this pre- 
cursor during one of its serendipitous spectral snapshots (Fig. 3), allowing us to get a good 
spectrum of it. The 3-250 keV spectrum (Fig. 4) is well fit by a single blackbody component, 
assuming N H = 6.69 x 10 22 cm" 2 , with kT = 10.4 ± 0.3 keV (xl = 1-06, 75 dof). By com- 
parison, both a simple power law model (xl = 4.02, 75 dof) and a thermal bremsstrahlung 
model (xl — 2.03, 75 dof) give unacceptable fits to the precursor. The spectrum shows 
no evidence for nonthermal emission in addition to the simple blackbody. Using this spec- 
tral snapshot, we can estimate the time-integrated blackbody fluence of this precursor to be 
(3.2 ±0.5) x 10~ 5 erg cm -2 , implying an energy of 3.8 x 1 41 d% n erg. Note, thi s spectral fit and 
fluence differ significantly from the results reported in iHurley et al. J2005h due to an error 
in our preliminary analysis of this precursor. (This error did not affect the rest of our pre- 
liminary analysis.) This current analysis resolves most of the apparent dis crepancy between 



our pr evious precursor fluence and that reported from ACS/INTEGRAL (IMereghetti et al. 
2005ah . 



4.2. Preflare 

Following the precursor, SGR 1806-20 appears relatively quiescent for 142 s until the 
main peak of the giant flare. While there is no obvious evidence in the lightcurve for emission 
during this period (Fig. 1), RHESSI snapshot spectra allow us to search for emission with 
better sensitivity than the lightcurves alone. We see no evidence of emission from SGR 1806- 
20 during this quiescent period. The 3a upper limit on the 3-250 keV flux is 1.4 x 10~ 7 
ergcm _2 s _1 , corresponding to a source flux below 1.7 x 10 39 d\ ergs -1 . 



4.3. Fast Peak 

There was somewhat less than one millisecond between the time that the first signs of 
the main peak became detectable above background by RHESSI and the time the instrument 
went into saturation. Fig. 5 shows the raw count rate at the start of the rise. We fit the early 
data where the instrument livetime was > 90% with an exponential growth curve, shown as 
well in Fig. 5. Normally, rear segment response at low energies on RHESSI is dominated by 
Compton reflection from Earth's atmosphere, but during the fast rise these photons would 
not have had a chance to reach the instrument, and therefore all segments are included in the 
lightcurve of Fig. 5. The best-fit e-folding time constant is 0.38 ± 0.04 ms. This rise time is 
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comp arable to the ~ 0.3 ms rise time reported for Swift-BAT before saturation (jPalmer et al. 
20051 ). but an order of magnitude fa ster than the 4.9 ms ri se time reported by particle detector 
observations that did not saturate ([Schwartz et al. 1120051 ) . However, we can also see in Fig. 5 
that RHESSI partly recovers during the t 2 6 = 638 — 639 ms period, suggesting that the onset 
of the main peak was not smooth, but had a significant drop in the incident rate ~ 2.5 ms 
after the st art of the rise. Th is drop is also evident in the Swift-BAT lightcurve at this 
same time (jPalmer et al.ll2005l ). One explanation of the apparent discrepancy between the 
rise time measurements is that RHESSI and Swift-BAT are characterizing the rise time of 
this initial fast peak, while the particle detectors are characterizing the main peak itself, 
following this initial fast peak. 

The fast rise time and short duration of this fast peak suggest that it represents a 
separate physical mechanism from the main peak itself. This idea is also supported by 
the spectrum of the fast peak. Fig. 6 shows the RHESSI count spectrum during the t 2 e = 
636.2 — 636.7 ms period. RHESSI did not get a direct snapshot spectrum during this period, 
so only a rough spectral analysis is possible. We can compare this count spectrum with 
the measured precursor spectrum with kT = 10.4 keV (solid line, Fig. 6). Here we plot 
the precursor count spectrum after the snapshot period (Fig. 4) for a direct comparison 
between the two count spectra outside snapshot periods. This fast peak spectrum appears 
harder than the precursor. However, this spectrum is much softer than the kT = 175 keV 
blackbody (dashed line, Fig. 6) that we measure for the main peak itself (next section). 
A kT no 20keV blackbody convolved through the rotation-averaged instrumental response 
gives a reasonable match to the measured count spectrum (dot-dashed line, Fig. 6). If we 
assume this component is a 20 keV blackbody, we can derive a rough ffuence for this 0.5 ms 
period during the fast rise of 6.6 x 10~ 7 erg cm -2 . 



4.4. Main Peak 

The Wind SST F & O count spectra of the main peak are shown in Fig. 7. As discussed 
in § 3, the Wind SST detectors measured the main peak spectral shape, but not the overall 
normalization. We concentrated our efforts on the F & O spectra during this analysis - 
the combined spectra of six detectors each. We exclude the so-called FT and OT coincident 
spectra due to the poorer statistics in this data mode, and uncertainties in the coincidence 
trigger criteria. In addition, the lowest bin in the F spectrum had to be excluded due to 
uncertainties concerning trigger threshold effects. In order to determine the best-fit spectral 
model, we developed a simulation mass model of the Wind 3D experiment, including a 
detailed detector and housing model, and a rough spacecraft model. Special care has been 
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taken to correctly model the passive material on all direct paths from the magnetar to 
the active detector materials and the immediate surrounds of the detectors to correctly 
account for absorptio ns and scatters. The mass model itself was developed in MEGAlib 



( IZoriauer et al. 1 120061 ) to allow both GEANT3 and GEANT4 Monte Carlo simulations in 
order to cross check results. Of specific concern was also the correct handling of electron 
tracks from Compton and pair interactions given the thin F & O detectors. However, since 
no calibration exists for this detection mode, some uncertainties remain. Given that most 
photons above roughly 50 keV interacting in one of the SST detectors will Compton scatter 
out of the instrument, the photon response matrix is strongly non-diagonal. As consequence, 
we did not attempt to create and invert a photon response matrix and produce model- 
independent spectra, but rather compared measured and simulated count spectra directly 
to determine the best match to the overall measured spectral shapes. 

We attempted to reproduce the observed count rate distributions with power law, ther- 
mal bremsstrahlung, and blackbody spectral models. We varied the input spectral parame- 
ter for each model (photon index, temperature) over a range of values, adjusting the overall 
normalization to best match the observed count spectrum. We verified that the range of 
spectral values bracketed the best-fit value for each spectral model. The best-fit power law 
and bremsstrahlung spectra were strongly rejected (xl — 4.2 and 6.9, 10 dof), and only the 
blackbody with kT = 175 ± 25 keV provided an acceptable fit (xl — 1-0, 10 dof). 

Given the best-fit spectral model from Wind SSTs, we can use the RHESSI PDs (Fig. 2) 
to constrain the overall fluence. The RHESSI PD was also modeled in GEANT3, assuming 
the best-fit blackbody model derived with the Wind SSTs above. The RHESSI PD data yield 
an overall fluence of (1.36 ± 0.35) erg cm -2 (for the blackbody integrated over all energies), 
implying an isotropic energy release in the main peak of 1.6 x 10 46 d\ Q erg. Given the RHESSI 
PD time resolution, the peak flux in the first 0.125s was 1 x 10 47 d\ ergs -1 . 

While the main peak is evident in only two RHESSI PD 0.125-s time bins, it shows clear 
indication of softening between these two intervals (Fig. 2). Assuming an input blackbody 
spectrum, we can use the two PD count rates to characterize the blackbody temperature 
independent of the Wind results. In the first and second time bins, the RHESSI PD charac- 
terizes the temperatures as kT = 230 ± 20 keV and kT = 170 ± 5 keV respectively. 



Ha j das et al. I (120051 ) performed a similar analysis on the IREM radiation environment 
monitor on the INTEGRAL satellite. IREM has three Si-diodes with crude spectral response, 
consisting of count rates read out every 60 s for events above six energy thresholds spanning 
85 keV to above 3 MeV. The IREM spectral analysis results in a cooling blackbody with 
kT = 230 ± 50 keV and an overall fluence of (0 .97 ± 0.50) erg cm -2 , consistent with our 



analysis above. By contrast, iMazets et al. I (120051 ) performed spectral analysis of the 7-ray 
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signal measured by the Russian spacecraft Coronas-F. This spacecraft was occulted by the 
Earth at the time of the main peak, but observed the event reflected off of the Moon. 
They simulated the response of scattering off the Moon, then folded this response with 
their detector response matrix. Their best-fit model for this spectral analysis is a powerlaw 
function (r = 0.7) with exponential cutoff (E a = 800 keV). However, from our own spectral 
analysis of the Wind SST data, we can rule out this spectral model with very high confidence 
(xl > 10, 11 dof). Therefore, given the quality of the Wind SST data, the consistency 
with the limited RHESSI PD data, and the consistent results with the IREM data, we are 
confident of the kT ~ 200 keV cooling blackbody spectrum of the main peak. 



4.5. Peak Decay 

When they came out of saturation within 1 s after the start of the main peak, the 
RHESSI spectroscopy detectors were measuring a peak count rate of ~280,000 cnt s -1 . Dur- 
ing the first few seconds, RHESSI recorded a dynamic and complex spectrum. What stands 
out most in the lightcurve (Fig. 8) is the pulsed tail, composed of both thermal blackbody 
emission and a nonthermal power law emission. Both of these spectral components are 
present immediately when RHESSI comes out of saturation, with the earliest snapshot at t26 
= 2.0 s. These two components, and their evolution throughout the pulsed tail, are discussed 
below. 

An additional component present as RHESSI emerged from saturation consists of strong 
emission extending up to 17 MeV, the upper limit of RHESSI's energy band. To our knowl- 
edge, this is the highest energy to which this or any other SGR flare has been observed. 
Fig. 9 shows the 0.4-10 MeV lightcurve for the first 300 seconds after the giant peak, includ- 
ing our fit to the background rate before and after the giant peak. RHESSI observes excess 
emission in the MeV range for ~60 s after the giant peak, which is better modeled as a power 
law decay than an exponential decay. Fitting the background-subtracted lightcurve with a 
function of the form oc t~ a yields a best-fit index a = 0.68 ± 0.04 (xl = 1.63, 39 dof). An 
MeV component was prev iously reported in t he pulse-averaged spectra from Konus-Wind 



observations to ~10MeV (iMazets et al.ll2005l ). though we do not see any significant MeV 
emission after t26 ~ 100 s. 

We verified that the high-energy component is not created by low-energy photons from 
the SGR arriving simultaneously with high-energy background photons and creating an arti- 
ficial "pileup" component at high energies. During the period of high-energy emission in the 
peak decay, the detector livetime averaged around 96% in the RHESSI rear segments. Using 
pileup-modeling software based on RHESSI ground calibrations and solar flare observations, 
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we find that the high-energy contribution of pileup should be about a factor of 10 lower than 
the high energy component observed. 

For the t26 = 1.71-22.06 s period, the 0.4-15 MeV spectrum (Fig. 10) can be fit by a 
power law oc E~~ r ) of photon index T = 1.43 ± 0.06, with an integrated fluence of 
(9.8 ± 0.1) x 10~ 5 erg cm -2 . For this fit, we excluded the band around the 0.511 MeV 
background line, which is difficult to model and subtract properly for this transient event. 
The power law model for this MeV component, shown in Fig. 10, is a good fit above 0.4 MeV 
(xl = 0.86, 30 dof). There is no sign of a turnover in this spectrum up to 17 MeV. Adding an 
exponential cutoff to the models marginally worsens the spectral fit, and results in a cutoff 
energy >50 MeV, well above the RHESSI energy range. 

Below 0.4 MeV, Fig. 10 shows a strong excess above our simple power law model. This 
excess is partly due to photons from SGR 1806-20 that scattered in the earth's atmosphere 
before reaching the spectroscopy detectors, and due in smaller part to a softer power law 
index at lower energies as revealed in the snapshot spectra (Sec. 4.6). While the snapshot 
spectra themselves show no sign of a change in the spectral index (~2.5) below 250 keV, the 
higher energy data in Fig. 10 show that the spectral index steepens (~1.5) at higher energies. 
The exact energy of this spectral break is only weakly constrained by the RHESSI data to 
be 0.5 ±0.2 MeV. 



4.6. Pulsed Tail 



After the main peak of the giant flare, RHESSI recorded a series of 51 pulsations with 
a period of 7.56 s (Fig. 1) , similar to the INTEGRAL, KONUS, an d Swift-BAT observations 
( IMereghetti et al.l l2005al ; iMazets et al.l 120051 ; iPalmer et al. 1 120051 ) . The pulse profile shows 
evidence for both spectral variations throughout the pulse, and evolution of the pulse shapes 
throughout the decay. The 20-100 keV pulse profiles show 3-4 peaks in their structure. 

In Table 1 we present the phase-integrated spectral evolution of the pulsed tail as seen 
through the snapshot spectra. Data were combined over the time periods presented, and fit 
in the 3-250 keV band (assuming Nh = 6.69 x 10 22 cm' 2 ) to determine the best spectral 
model. For these pulse-integrated spectra, the best-fit model is a two-component blackbody 
plus power law (photon index T). Thermal bremsstrahlung and power law models give 
unacceptable fits, both alone and combined as a two-component model. A blackbody model 
alone gives marginal fits, improved significantly with the addition of the nonthermal powerlaw 
component. We do not see any significant sign of an exponential cutoff in this power law 
component below 250 keV, the upper end of our snapshot data. After t 2 6 ~ 246 s, the 
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power law photon index is not strongly constrained in the individual 40-s spectra though 
the component is still significant; therefore, we have fixed the index at its best-fit value for 
this time period, T = 2.1. An example spectral fit is shown in Fig. 11, the first snapshot 
spectrum after RHESSI emerged from saturation. 

In our preliminary analysis we could not conclusively distinguish between the therma l 



bremsstrahlung and blackbody models, which both gave marginal fits (IHurley et al. 1 120051 ) . 
The addition of the absorption column and the power law component improved the blackbody 
fits to the point of strongly distinguishing the models. 

The blackbody component appears to be present in the tail emission immediately after 
RHESSI comes out of saturation, with an initial temperature kT = 11.5 keV, which drops 
steadily during the evolution of the pulsed tail (Table 1). The total integrated fluence of the 
blackbody component of the pulsed tail is (2.6 ± 0.2) x 10~ 3 erg cm -2 . 

The power law component is also present immediately after RHESSI comes out of sat- 
uration, and lasts throughout the pulsed tail phase. The photon index, initially T = 1.71, 
appears to soften to V = 2.7 then harden to T = 2.1 over the evolution of the tail, but it 
is not clear that this evolution is strongly significant. The total integrated fluence of the 
nonthermal component for the pulse tail phase, in the 3-100 keV band, is (2.9 ± 0.5) x 10~ 3 
erg cm -2 . 

Combining these two components, the total fluence measured from the pulsed tail is 
(5.5 ± 0.6) x 10~ 3 erg cm -2 , implying an energy release of 6.7 x 10 43 d\ Q erg, roughly equally 
divided between thermal and nonthermal emission. 

Fig. 12 shows the average 7.56-s pulse shape, 20-100 keV, integrated over the pulsed 
tail. The pulse profile is dominated by 3-4 separate peaks, with an overall large pulse 
fraction. Fig. 12 also shows the phase-resolved best-fit blackbody temperature, which varies 
throughout the pulse. During the phase-resolved profile, the ratio of 3-100 keV power law 
flux to the total blackbody flux stays flat, with the energy emission nearly equally divided 
between the two components. 



4.7. Afterglow 

ACS/INTEGRAL observations of the SGR 1806-20 giant flare showed excess counts in 
the period following the pulsing phase, t 2 6 ~ 400-4000 s, peaking around t26 ~ 600-800 s. This 



excess was interpreted as afterglow hard X-ray emission from the SGR (jMereghetti et al. 



2005al ). The ACS observations suggest that the integrated fluence in this afterglow emission 
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is comparable to the integrated emission in the pulsing tail itself. RHESSI had a direct view 
of SGR 1806-20 for 600 s following the pulsing tail before the satellite moved behind Earth's 
shadow. We have searched the RHESSI snapshot spectra during this period for evidence 
of afterglow emission. For t 26 ~ 400-1000 s, we can set a 3a upper limit on the 3-200 keV 
ffuence from SGR 1806-20 of 1.4 x 10~ 4 ergcm" 2 . This upper limit is a factor of 60 below 
the ffuence expected based on the ACS observations. Given our snapshot spectra, we are 
effectively chopping between source and background. Therefore, our analysis is not sensitive 
to potential background variations as is the ACS lightcurve analysis. Activation of the BGO 
crystals comprising the INTEGRAL ACS detectors was considered as a possible explanation 
for the excess counts, but it is not clear that this would be consistent with the observed 
lightcurve of the ACS afterglow (Mereghetti, private communication). 



Discussion 



The precursor may hold some tantalizing clues to the origin of this giant flare. The 
unusual nature of this burst, and its occurence soon before the main peak, suggest a direct 
connection between the two events. Ind eed, a precursor was also present for the giant flare 
from SGR 1900+14 (jHurley et al.lll999l ). The average luminosity during the precursor was 



4 x 10 41 d\ Q ergs -1 , corresponding to ~2000L £ ; (M , which is typical for SGR bursts. In 
Fig. 13 we show the phase of this precursor relative to the subsequent pulsed tail - the 
precursor occurs during one of the lowest phases of the pulse profile. Our measured rise 
time of this precursor, 27ms, is much longer than the ~0.03ms timescale expected for 
magnetospheric realign ment, but consistent with t he expected 10 ms timescale for crustal 
slipping after cracking ([Thompson fc Duncan 1 120011 ). This timescale is strong evidence that 
the precursor originate d from a fracture prop agating in the crust of the neutron star, just as 
in the main peak itself (jSchwartz et al.ll2005l ). The total energy of this precursor, 3.8 x 10 41 
d\ Q erg, is comparable to the maxi mum elastic potential energy the magnetar can store in 
its crust before cracking, ~10 42 erg (IThompson fc Duncan 1 1200 ll ). which might suggest that 
the precursor corresponds to a global cracking and realignment of the crust. However, the 
relatively long 1-s duration of the precursor (5x the main peak of the giant flare itself) 
and the multi-peaked lightcurve suggests instead that the precursor is being powered by 
repeated injections of energy by re alignment of the magnetic fie ld in the core, with a typical 
times cale on the order of 200 ms (IThompson fc Duncan 1 1200 ll ). similar to the main peak 
itself ( jSchwartz et al. 1120051 ) . Given the energetics, this scenario suggests that the precursor 
corresponds to an energetically small crustal fracture, followed by repeated energy injections 
from a relatively small realignment of the core. This scenario is supported by two more pieces 
of evidence. First, our precursor spectrum is purely blackbody. If there were significant 
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twist ing of the magnetosphere we would have expected a nonthermal tail ([Thompson et al. 
20021 ). Second, the drop in flux by over a factor of 200 after the precursor (during the preflare 
period) also suggests that the magnetosphere was not significantly twisted by the precursor 
event itself. 

We measured the peak flux in the first 0.125 s of the main peak to be 1 x 10 47 d\ Q 
ergs -1 , an astounding l$h Edd . In our previous paper we showed that this luminosity and 
the measured temperature of kT ~ 200 keV are consistent for blackbody emission from a 
spherical surface of radius R ~ 10 km (see below). We have measured the isotropic energy 
release in the main peak of 1.6 x 10 46 df erg. This is an awesome amount of energy for this 
source. For comparison, given the 7.56 s period the rotational kinetic energy of SGR1806- 
20 (the power supply for radio pulsars) is on the order of E spin = \l£l 2 = 3 x 10 44 erg. 
Given that the energy release in the main peak is two orders of magnitude greater than this 
rotational energy of the star, it is even more amazing that t here was no measurable jump in 
the spin frequency after the giant flare (IWoods et al. 1 120061 ). 



The energy release of the main peak is comparable to the maximum energy that could be 
stored in a twisted magnetosphere, thus energetically t his giant flare would be consistent with 
global magnetospheric untwisting (jliurley et al. 2005 ). However, in this scenario one would 
expect a growing spin-down rate before the giant flare, whereas it was actually decreasing 
in the months prior, and a significant drop in the spin-down rate after the giant flare, likely 
large r than the observed decreased which appears consistent with the trend before the giant 
flare (IWoods et al. 1120061 ). Thus we are left with the conclusion that this giant flare represents 
a large-sc ale crustal instability in the star, driven by the unwinding of the toroidal field inside 
the core ( Thompson fc Duncan|[2~00ll ). The main peak taps only a fraction of the 10 48 - 10 49 
erg of magnetic energy stored in the core of the SGR. 

This conclusion appears consistent with the strong nonthermal emission during the 
peak decay and throu ghout the pulsed tail, w hich indicates large magnetospheric twisting 
during the main peak (IThompson et al. 1 120021 ) . The nonthermal emission during the pulsed 
tail, with spectral indices ~ 2.1-2.7 and extending >250keV, is likely easily explained by 
electron c yclotron scattering with in an extended corona, which can reach photon energies 
>100keV (IThompson et al. 1 120021 ). However, the second nonthermal component during the 
peak decay, with a spectral index of 1.43 and extending >17MeV with no sign of a spectral 
cutoff, is more difficult to explain. It clearly must derive from a separate mechanism than 
the electron cyclotron scattering. Ion cyclotr on scattering is only expected to extend into 
the tens of keV range (IThompson et al. 1 120021 ). We speculate that this component perhaps 
arises from a highly extended corona, driven by the hyper-Eddin gton luminosities, where 
synchrotron emission is no longer efficient at cooling the electrons (IFeroci et al.ll200ll ). This 
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scenario seems consistent with the lack of a clear pulsation from the SGR for over a full 
rotation period, nearly 10 s, following the main peak. 

Fig. 13 shows that there is significant evolution of the pulse profile over the course of the 
pulsed tail. Each panel shows the average pulse profile integrated over ten consecutive pulses. 
Especially significant is the increase of the peak at phase ~0.65 through the evolution of the 
tail. Evolution of the pulse profile demonstrates con tinuing magnetospheric realignment 
during the course of the pulsed tail. iFeroci et al. I (120011 ) have studied a similar phenomenon 
in the pulse profile of SGR 1900+14 following a giant flare. 

The phase-averaged luminosity during the 400 s pulsed tail corresponds to 300-2000 Lsdd- 
In our previous paper, we showed that the pulse- averaged 20 -100 keV flux light curv e of th e 
pulsed tail is well modeled by the trapped fireball model of iThompson &: Duncan I (120011 ). 
where flux oc ( 1 — ^fz^ 1 ^".^ with an evaporation time t evap = 382 ± 3 s, and index a 



0.606 ± 0.003 (IHurley et al. 1 120051 ). This index is physically signific ant, being close to the 



expec ted value a = f for a homogeneous, spherical trapped fireball (IThompson fc Duncan 



200ll ). Within this trapped fireball model, we can derive a rough bound on the magnetic 
field by requiring the magnetic field to be strong enough to c onfine energy radiated by th e 
trapped fireball, B dipole > 1.7 x lQ^fJ^ir + ,/2] 3 dm G ( IThompson fc Duncanlll995h . 



We observed blackbody components, of different temperatures and luminosities, from 
the precursor, main peak, peak decay, and pulsating tail. For a source distance d, and 
surface gravitational redshift z, we can convert measured luminosities and temperatures into 
effective blackbody radii of the emission regions, R = {-^)*{jt^)- I n Fig- 14 we present our 
derived blackbody radii for the various stages of the giant flare, assuming d = lOkpc, and 
z = 0.30 for a neutron star (4r ~ 1 '^ Is ° 1 )- The first remarkable feature about these radii 

V it lOfcm ' 

is that they all roughly agree throughout the stages of this giant flare, from the precursor 
(14 km) and the main peak (18 km), through the main peak decay (17 km) and the average 
throughout the pulsing tail (11 km). These agree remarkably well given the variation in 
temperature and luminosity throughout. In addition, we can turn this around to point out 
that, in terms of the uncertain distance to SGR 1806-20, d < 10 kpc would be more consistent 
with a canonical neutron star radius of R ~ 10 km than d = 15 kpc if the thermal emission 
is originating from the stellar surface as opposed to an extended corona. 

The fast peak just prior to the main peak remains a bit of a puzzle. The average lumi- 
nosity during the 0.5 ms before RHESSI saturates was 1.6 x 10 43 d\ Q ergs -1 , corresponding to 
~ l^PljEdd- The peak in the lightcurve and the relatively soft spectrum clearly distinguish 
this as separate from the main peak itself. However, the rise time (0.4 ms) and duration 
(2.5 ms) of this fast peak are too fast for global crustal slippage or realignment of the core 
field (IThompson Duncan 1 120011 ). These times could imply a very localized crack in the 
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crust (~ 0.4km), but this fast peak is much shorter than typical SGR flares (and the pre- 
cursor), which are attributed to the same physical mechanism. Furthermore, if we assume 
the spectrum really corresponds to a blackbody of kT ~ 20 keV, then we can also derive a 
blackbody radius of roughly 21km for this stage of the flare, i.e. a global event, presumably 
inconsistent with a localized crustal crack. If we divide this radius by the rise time we can 
derive a lower limit on the thermal diffusion speed, corresponding to ~0.2c. This fast peak 
appears to only be consistent in timescale and sizescale with a realignment of the magneto- 
sphere immediately before the main peak itself. While the energy of this fast peak is small 
compared to the main peak, its fast timescale and its occurrence milliseconds prior to the 
main peak suggest it plays some critical role in the giant flare. 

The giant flare of SGR 1806-20 represents one of the most outstanding events in X-ray 
and 7-ray astronomy over the past three decades, even when compared to the giant flares 
of SGR 0525-66 and SGR 1900+14. Even though this giant flare presents an extreme energy 
output from the SGR, the basic energetics and timescales involved are well understood in 
term of the magnetar model. Indeed, given the spin energy of the star and the extraordinary 
super-Eddington luminosities of the giant flare, SGR 1806-20 has once more presented very 
strong evidence in favor of the magnetar model. 

The authors are grateful to G. Hurford, H. Hudson, and S. Krucker for useful discussions. 
KH is grateful for support under the NASA Long Term Space Astrophysics program grant 
NAG5-13080. 
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Fig. 1. — The RHESSI 20-100 keV background-subtracted time history of the giant flare, 
plotted with 0.5-s resolution. The main peak begins at 0.64 s, where the RHESSI detectors 
are saturated and effectively dead. RHESSI recovered ~0.5s later to observe the rest of 
the giant flare in detail. In this paper we analyze six separate stages of this flare: (i) 
precursor, (ii) preflare, (iii) fast peak, (iv) main peak, (v) peak decay, (vi) pulsed tail, and 
(vii) afterglow. In this energy range, the time history is modulated by RHESSI's 4-s spin 
period; otherwise, there are no other long-term induced modulations (instrument repointing, 
saturation, etc.). Therefore, other than the main peak this represents the true lightcurve of 
the giant flare event. 
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Fig. 2. — The RHESSI particle detector was able to measure the incident flux with 0.125-s 
time resolution in two energy channels determined by thresholds in the electronics: >50 keV, 
and >620 keV. These data indicate significant emission above 620 keV for ~0.25 s, during 
the main peak, and softening of the spectrum during its evolution. 
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Fig. 3. — Precursor 20-100 keV lightcurve, in ^ s time bins. The precursor has an exponen- 
tial rise time of 27ms, and a fall time of 110ms (solid lines). The period of the precursor 
snapshot spectrum (Fig. 4) is shown by the dashed lines. 
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Fig. 4. — RHESSI snapshot spectrum of the precursor, 3-150 keV. This spectrum is 
background-subtracted, and the solid line shows the best-fit absorbed blackbody spectrum 
convolved through the RHESSI response matrix. 
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Fig. 5. — Total RHESSI count rate during the rising edge of the giant flare. Beyond the 
dashed line, instrument deadtime becomes significant; in fact, the peak in count rate at 638.5 
ms represents a dip, not a peak, in the true flux from a partial recovery in livetime. The 
instrument is almost completely paralyzed during the periods of low count rate after 637 ms. 
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Fig. 6. — RHESSI count spectrum during the fast peak just prior to the main peak. RHESSI 
did not get a direct snapshot spectrum during this period, so only a rough spectral analysis 
is possible. The measured precursor spectrum (outside of the snapshot period) is shown for 
comparison (solid line). Also shown for comparison are two blackbody spectra convolved 
with the instrument response matrix: kT = 175 keV corresponding to the main peak tem- 
perature (dashed line), and kT = 20keV (dot-dashed line). The normalizations on the three 
comparison spectra are arbitrary. 
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Fig. 7. — Count spectra of the main peak of the Wind SST F and detectors. The crosses 
show the measured spectra and the lines the best-fitting simulated spectra, a blackbody with 
the temperature 175±25 keV (xl = 1.0, 10 dof). 
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Fig. 8.— The RHESSI 20-100 keV time history in 0.125 s bins as RHESSI comes out of 
saturation, measuring the main peak decay and the transition to the pulsed tail. 
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Fig. 9. — The RHESSI 0.4-10 MeV lightcurve, 4.07s bins (spacecraft rotation period), show- 
ing the high energy excess in the decay phase as RHESSI comes out of saturation. The blue 
curve shows the best fit to the underlying background rate. The red curve shows the best-fit 
power law decay to the MeV excess rate. The green lines mark the time interval used for 
spectral analysis (Fig. 10). 
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Fig. 10. — RHESSI MeV spectrum of the high energy excess, integrated over the 20-s interval 
shown in Fig. 9. The blue curve shows the best-fit power law model of the 0.4-17 MeV spec- 
trum (red points). The spectral bin including the 0.511 MeV background line was excluded 
from this fit. 
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Fig. 11. — The first RHESSI snapshot spectrum after coming out of saturation, at t26 = 
2.0s. The solid line shows the best-fit absorbed blackbody + powerlaw model convolved 
through the RHESSI response matrix. The absorbed blackbody (dashed line) and absorbed 
powerlaw (dotted line) are also shown separately. 
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Fig. 12. — Average pulse shape, 20-100 keV (top), integrated over the entire 400s pulsed 
tail. The bottom panel shows the phase-resolved blackbody temperature. 
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Fig. 13. — Evolution of the 20-100 keV pulse profile during the pulsed tail. Each frame shows 
the pulse profile integrated over 10 successive pulses, from the start of the pulsed tail (top) 
to the end (bottom). In the top panel, we also show the phase of the precursor relative to 
the pulsed tail. 



31 




Fig. 14. — The effective blackbody radius for various stages of the giant flare, assuming a 
distance of lOkpc and a gravitational redshift of z — 0.30. 
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Table 1: Pulsed tail integrated blackbody + power law spectral fits (3-250 keV). 



t26 


kT 6fe 


BB Fluence 1 


r 


PL Fluence 2 


Xl (dof) 




[keV] 


[10~ 4 ergcmT 2 ] 




[10~ 4 erg cm -2 ] 




2-6 s 


11.0 ±0.2 


3.1 ±0.4 


1.71 ±0.09 


0.6 ±0.2 


1.11 (73) 


6-46 s 


9.7 ±0.3 


5.0 ±0.9 


2.34 ±0.09 


6.3 ±2.0 


1.23(73) 


46-86 s 


9.6 ±0.4 


3.2 ±0.7 


2.58 ±0.14 


5.4 ±2.2 


1.29(73) 


86-126 s 


9.1 ±0.3 


3.5 ±0.7 


2.49 ±0.16 


3.8 ±2.0 


0.98 (73) 


126-166 s 


9.4 ±0.3 


3.4 ±0.7 


2.64 ±0.18 


4.4 ±2.4 


1.01 (73) 


166-206 s 


8.4 ±0.4 


2.4 ±0.7 


2.70 ±0.15 


5.0 ±2.5 


1.20(73) 


206-246 s 


7.7 ±0.5 


1.5 ±0.6 


2.50 ±0.20 


2.3 ± 1.7 


1.06(73) 


246-286 s 


5.9 ±0.2 


1.3 ±0.3 


2.1(fixed) 


0.7 ±0.1 


1.53 (74) 


286-326 s 


5.5 ±0.2 


1.1 ±0.2 


2.1(fixed) 


0.6 ±0.1 


1.16(74) 


326-366 s 


4.4 ±0.2 


0.9 ±0.2 


2.1( fixed) 


0.4 ±0.1 


0.98 (74) 


366-406 s 


3.5 ±0.4 


0.3 ±0.1 


21(fixed) 


0.06 ±0.07 


0.86 (74) 



1. Blackbody fluence integrated over all energies. 

2. 3-100 keV fluence. 



